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1. Introduction 

Carbon nanotubes [1, 2] are one of the most promising materials for future 
nanoelectronics due to their unique electrical transport properties [3] . In perfect single- 
walled carbon nanotubes (SWCNTs), baUistic electron conduction has been observed [4] 
provided inelastic processes can be neglected. In this coherent regime, however, the 
presence of defects, dopants or other impurities will particularly affect the performance 
of the device [5], as it will determine the transport in nanotubes from a ballistic regime 
to either weak or even strong localization regimes, where the resistance will increase 
exponentially with the length of the system. 

To fully achieve control on the properties of carbon nanotubes, it is thus essential 
to perform accurate analyses of realistic disorder on carbon nanotubes-based devices. 
Although simulation tools based on Density Functional Theory (DFT) provide great 
accuracy and are particularly suitable for the study of defects in a wide range of 
materials, their applicability is limited by their poor scaling, which leads to a huge 
computational effort when the number of atoms in the system is larger than a few 
hundred. This prevents from using standard ab initio techniques to analyze directly the 
properties of electronic devices in the mesoscopic scale. 

To overcome this limitation, we have combined the use of standard (no linear- 
scaling) first principles methods to extract the electronic properties of single defects 
in carbon nanotubes with the use of 0(N) Green's functions techniques to study the 
transport properties of long tubes with a random distribution of those defects. This 
approach allows to calculate characteristic transport lengths such as the localization 
length for different densities of defects, providing an interesting tool to tune the 
transport properties of carbon nanotubes, for instance, by a controlled creation of the 
defects [6] through irradiation, and was able to predict the existence of the Anderson 
localization regime in these carbon systems at room temperature as long as the electronic 
phase coherence is preserved [6, 7]. The emergence of the Anderson localization 
phenomenon in chemically doped carbon nanotubes has been theoretically predicted 
in other works [8,9], in which a similar approach based on a tight-binding modeling 
from ab initio calculations of single B or N impurities was used to characterize the 
transition between different transport regimes. 

2. Theoretical framework 

Our analysis of the transport properties of the defected carbon nanotubes was separated 
into two steps. First, we studied the perturbation that single, isolated defects cause -at 
zero temperature- on the geometry and the electronic properties of an otherwise ideal 
tube. In a further step, a study of the transport properties of nanotubes with a more 
realistic disorder has been performed, calculating the conductance along the nanotubes, 
both at zero and finite temperature, for a statistically meaningful sample of nanotubes 
with different concentrations of randomly distributed defects. 
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Our approach combines first-principles calculations to obtain the effective one- 
electron Hamiltonian for the defected nanotubes with non-equilibrium Green's functions 
techniques, which allow the calculation of the conductance for very long -up to microns- 
tubes. 

2.1. Geometry optimization 

In our calculations we have used the ab initio Firebair96 code [10 12], that provides fast 
and accurate molecular dynamics calculations of the geometries and electronic properties 
for a wide range of systems. This method uses localized orbit als, generated by solving 
the atomic problem within the DFT-LDA and the pseudopotential approximations. 
Norm-conserving pseudopotentials are used, as well as the Ceperly-Adler form of the 
exchange-correlation potential as parametrized by Perdew and Zunger [13]. The forces 
on each atom are calculated by means of a variation of the Hellmann Feynman theorem, 
and molecular dynamics are used to obtain the lowest-energy atomic configuration. A 
self-consistent implementation of the Harris functional allows to avoid the much time- 
consuming calculation of the electronic charge density p(r), which is substituted by the 
calculation of the orbital occupancies n^^ of the 'fireball' orbitals 0^, defined by 

P(r-) = E^mI0>")I' • (1) 

This approach allows us to optimize the electronic charge density p according to the 
chemical environment of the atoms. The different values of are then obtained by 
imposing a self-consistent condition on the orbital occupancies. 

In the simulations with Firebair96 an sp^-basis set of orbitals was used with a 
cut-off radius of 2.15 A. The supercell contained 12, 10 and 6 unit cells for the (5,5), 
(7, 7) and (10, 10) nanotubes, respectively (the unit cell containing 20, 28 or 40 atoms 
for each tube) . Four special t-points uniformly distributed along the first Brillouin zone 
of the nanotubes have been included. Convergence of the calculations is achieved for 
forces acting over each atom being less than 0.01 eV/A and changes on the total energy 
of the system are below 10~^ eV; the time step was 0.5 fs. 

2. 2. Transport formalism 

The calculation of the conductance has been performed within the Landauer- 
Biittiker formalism, more concretely in the Keldysh (non-equilibrium Green's functions) 
formalism, which is particularly convenient when using a localized basis set. The 
differential conductance g between two systems in the limit of zero temperature and 
low voltages is given by [14, 15] 

9{E) = ^Tr[Dt,{E)f,,pf^{E)f,,D^,{E)pf^{E)] (2) 

where Tr[0] represents the trace of operator O; T12 describes the coupling between 
systems 1 and 2; Pn{E) and P22{E) are the density of states matrices associated with 
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the decoupled (T12 = 0) systems 1 and 2. Multiple scattering effects are included by 
means of the retarded and advanced denominator functions D\ 

J^n — [-^ ~ ^12^22 ^21^11 J 

= [I-f,,GTfi2GT]-' ■ (3) 

We are interested in calculating the conductance of a defected nanotube, but 
we are not including in this analysis the effect of realistic contacts. The simulation 
geometry will thus consist of a device region where the nanotube containing the defects 
is connected to two semi-infinite perfect tubes (L and R) that act as left and right 
electrodes (see Figure [1]); the density of states of the electrodes is calculated using 
standard decimation techniques. The conductance along the defected region of the 
nanotube is calculated by an iterative procedure, starting on the left electrode L, and 
including at each step a random number of ideal layers of nanotube (shown as . . ' in 
Figure ^ as well as the block (supercell) of the nanotube containing the defect, and 
whose most stable reconstruction has been obtained using Firebair96 (labeled as Ci and 
C2 in the same figure). After the inclusion of a new defect, the conductance between 




L q R 

Figure 1. (color online) Nanotube scheme as considered to perform the conductance 
calculation: L and R are the semi-infinite ideal nanotubes that act as left and right 
electrodes; the central region is the defected nanotube, consisting of one or more blocks 
of ideal ('...') or defected (ci, C2, ...) nanotube. 

the 'left side' of the system (consisting of the left lead and the increasing number of 
ideal and defected layers of nanotube) and the right electrode (which remains invariable 
during the process) is calculated. In this way, we can calculate the conductance as 
a function of the number of defects included (and thus of the length of the tube) for 
different random configuration of defects distributed along the nanotube with a given 
mean distance d between them. The length L of the tubes is obtained from L = N - 
where is the number of defects in the nanotube. 

The main advantage of this linear-scaling numerical procedure is the fact that the 
size of the matrices we need to invert at each step does not increase with the size of 
the system, but remains constant and equal to the size of the defected region of the 
nanotube. In this way, we can calculate the conductance for mesoscopic-size tubes with 
different concentrations of defects without a very expensive computational effort. 

The hopping or coupling between the layers of nanotube is assumed to be the 
coupling between the bulk layers of ideal nanotubes. For this reason, it is essential that 
the number of layers included in the supercell of the ab initio simulations is large enough 
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so that the surface layers of the superceU are not being affected by the perturbation 
caused by the defect and can be regarded as 'ideal' bulk layers. 

At zero temperature, the conductance for an infinitesimal voltage is calculated just 
by the evaluation of the differential conductance g{E) at the Fermi level Ef. However, 
at finite T there will be more states accessible to the electrons and the calculation of 
the conductance will be performed by means of the following expression [16]: 

G^J^jy-^),<E)iE, (4) 

where 

ME) = (5) 

is the Fermi distribution function. In this way, the total conductance at a given 
temperature is calculated as the total contribution from all the electronic states included 
in the thermal energy window opened by the finite temperature. 

3. Analysis of a single defect 

Theoretical simulations [17, 18] suggest that Ar^ irradiation on carbon nanotubes will 
mainly create monovacancies and divacancies along the nanotube. We have therefore 
studied the monovacancy and the two possible orientations of the divacancy, which we 
will refer to as vertical (if the orientation of the chain of missing atoms is perpendicular 
to the nanotube direction) and lateral (if the orientation is partially parallel to the tube 
axis) . 

3.1. Monovacancies 

The structural and electronic properties of vacancies-related defects in CNTs and other 
graphene-based materials have been extensively analyzed in the literature (see, for 
instance, [19-24]). The extraction of one single atom in the nanotube leads to a 
metastable state, that propitiates recombination of nearest neighbors to fill the empty 
space left by the vacancy. 

3.1.1. Effect on the geometry After the relaxation, we found that new bonds are formed 
between the neighbors of the extracted atom, favoring three-atoms against two-atoms 
coordination. This well-known recombination of atoms is found in the three types of 
armchair CNTs analyzed in this work. 

The tube diameter is also affected by the extraction of the single atom. It decreases 
in the zone of the defect and reaches its maximum value at a distance of 5 A from the 
defect, recovering the diameter of the ideal tube at the surface layer of the supercells 
used in our simulations. Table [H shows diameter (d) for ideal tubes as well as minimum 
and maximum diameter for defected ones. 
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Figure 2. (color online) Relaxed geometries of the (5,5) (top left), (7,7) (top right) 
and (10,10) (bottom) SWCNTs with a single monovacancy. Only the area around the 
defect (6 layers) is shown. 

Table 1. Maximum and minimum diameters for (5,5), (7,7) and (10,10) nanotubes 
with a single monovacancy and comparison of AEmonovac foi" the relaxed geometries. 



CNT 


dideal (^) 








(5,5) 


6.97 


6.56 


7.34 


16.62 


(7,7) 


9.76 


9.36 


10.20 


17.41 


(10,10) 


13.89 


13.86 


14.76 


17.42 



3.1.2. Energetics The analysis of the energetics of the monovacancy has been 
performed for each CNT by means of the following difference of total energies, that 
allow us to compare total energies for systems containing a different number of atoms: 

^-^monovac = | -^ideal {-^monovac ~l~ -^atom) \ •> (6) 

where Ei^eai is the total energy for an ideal tube calculated for a supercell containing 
the same number of layers than the supercell for the defected tube, Ej^onovac is the total 
energy of the relaxed nanotube containing the monovacancy, and Eatom is the energy 
for the free atom. AEj^ionovac for each tube is also shown in Table [H 

3.2. Divacancies 

Divacancies are created extracting a second atom form the nanotube, leading to two 
possible orientations of the divacancy depending on which atom is extracted: removing 
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Figure 3. (color online) Relaxed geometries of the (5,5) (top left), (7,7) (top right) 
and (10,10) (bottom) SWCNTs with a vertical divacancy. Only the area around the 
defect (6 layers) is shown. 

the atom placed just above (or below) the first atom extracted will create a vertical 
divacancy (Figure [3]), whereas the lateral divacancy is formed by the extraction of the 
atom located to one side of the monovacancy (Figured]). 

We have found a similar reconstructed geometry for the three CNTs analyzed: 
atomic relaxation leads to a new geometry formed by two pentagons and an octagon. 
As for the monovacancy case, the diameter decreases locally at the zone containing the 
defect, but the effect is much more pronounced for the case of the lateral divacancy (see 
Table El). 

3.2.1. Effect on the geometry 

3.2.2. Energetics The energy formation of the divacancies has been calculated in a 
similar way to the one performed in the case of the monovacancy. Now, 

^E divac — I E Ideal {Edivac H" 2 X E^iq^ji^ | . (7) 

According to the calculated energy formation, lateral divacancies are about 2 eV more 
stable than vertical ones for all the analyzed nanotubes. 

We can also estimate the energy for an ideal system with two monovacancies 
infinitely separated using 

E2vac = E ideal — "^{E ideal — Eyac) • (8) 
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Figure 4. (color online) Relaxed geometries of the (5,5) (top left), (7,7) (top right) 
and (10,10) (bottom) SWCNTs with a lateral divacancy. Only the area around the 
defect (6 layers) is shown. 

Table 2. Maximum and minimum diameters for (5,5), (7,7) and (10,10) nanotubes 
with a vertical and lateral divacancy and comparison of A for the relaxed geometries. 



CNT 


dideal (^) 


dmin (A) 


d^a. (A) 


A (eV) 


(5, 5) vertical divacancy 


6.97 


6.99 


7.19 


4.3 


(5,5) lateral divacancy 


6.97 


6.64 


7.28 


6.5 


(7, 7) vertical divacancy 


9.76 


9.72 


9.99 


4.8 


(7, 7) lateral divacancy 


9.76 


9.36 


10.02 


6.7 


(10, 10) vertical divacancy 


13.89 


13.88 


14.07 


4.2 


(10, 10) lateral divacancy 


13.89 


13.72 


14.65 


6.7 



Comparing this energy difference AEdivac with the previously calculated AEmonovac 

^ — I divac 2 X AEjYiQjiQydQ^ , (9) 

we have found that formation of a lateral divacancy is between 6 and 7 eV more stable, 
for all CNTs studied, than formation of two monovacancies. Energy formation of vertical 
divacancies is for all tubes between 4 and 5 eV smaller than that of the lateral one when 
compared to the case of two isolated monovacancies. Thus, lateral divacancies formation 
on CNTs is preferred over that of vertical ones, and especially than formation of isolated 
monovacancies. Results are summarized in Table [2l 
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3.3. Scattering of a single defect at zero temperature 

Although a large number of theoretical works have been performed to analyze the 
scattering produced by single defects (either topological or substitutional impurities), 
it is worth mentioning that a full relaxation of the atomic region around the defect 
is particularly important in the case of vacancies, where the energetically more stable 
reconstructed vacancy will present a different electronic structure than the ideal (non- 
reconstructed) one, leading to a very different behavior of the conductance and other 
transport properties, especially for energies close to the Fermi level. 

In the case of a tt — tt tight-binding model [25], a monovacancy in an armchair 
CNT will produce a drop of almost one conductance quantum. Go = 2e^//i, right at 
the Fermi level, with the electron-hole symmetry fully preserved. The conductance of 
a single ideal monovacancy, that is, in the metastable non-reconstructed state, will as 
well lead to a drop of almost 2e^//i, not at the Fermi level but somewhat below at the 
valence band, as shown in Figure [5] (where the Fermi energy has been set to zero). This 
result agrees with previous ab initio calculations [26,27]. 




_l I I I L 

-1 1 

Energy (eV) 



Figure 5. Conductance as a function of energy for the (10,10) SWCNT with a single 
monovacancy, before relaxation. 

However, the situation is very different in the case of the most stable reconstruction. 
Figure [6] shows the conductance as a function of the energy for the (5,5), (7,7) 
and (10, 10) nanotubes with a single monovacancy after relaxation. Rearrangement 
of atoms has propitiated the saturation of dangling bonds, and the reconstructed 
nanotube presents a conductance much more similar to the ideal case for low energy 
than that of the metastable monovacancy. Particularly, there is a very small drop in 
the conductance for energies in the vicinity of the Fermi level, and the narrow peaks 
in the conduction band due to the unsaturated a bonds that are present for the non- 
reconstructed geometry have disappeared. As expected, the drop in the conductance 
decreases with increasing diameter of the nanotubes, being practically negligible for the 
(10, 10) tube. This is in good agreement with previous calculations [27,28] 

The drop of the conductance is much larger for both orientations of the divacancy. 
In all cases a maximum drop of about one conductance quantum is found for the two 
orientations of the divacancy. This is also in agreement with previous works [23] . 
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Figure 6. Conductance as a function of energy for the (5,5) (top left), (7,7) (top 
right) and (10,10) (bottom) SWCNTs with a single monovacancy after relaxation. 
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Figure 7. Conductance as a function of energy for the (5,5) (top left), (7,7) (top 
right) and (10,10) (bottom) SWCNTs with a single vertical divacancy after relaxation. 



As lateral divacancies are much more stable than the vertical ones, and due to small 
drop in the conductance caused by monovacancies, the results we present include only 
lateral divacancies. 
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Figure 8. Conductance as a function of energy for the (5,5) (top left), (7,7) (top 
right) and (10,10) (bottom) SWCNTs with a single lateral divacancy after relaxation. 



4. Statistical analysis 

Once the transport properties of the isolated defects in the infinite nanotube have been 
characterized, we are interested in calculating the conductance along a more realistic 
nanotube in which defects can be randomly distributed along and around the axis of 
the tube. Lateral divacancies are distributed along the tube with the distance between 
consecutive defects presenting a uniform random distribution between and 2d -d thus 
being the mean distance between defects for a particular defect density. The calculation 
of the resistance (as the inverse of the differential conductance) for each particular 
random distribution of defects has been performed as explained in Section 12.21 

For each particular distribution of defects we can calculate the conductance as a 
function of energy according to In Figure [9] we plot the conductance around a small 
energy region around the Fermi level for the (5, 5) nanotube and a mean distance between 
defects of d = 10.1 nm for different number of defects included in the nanotube. The 
nanotube conductance presents strong fiuctuations, with a maximum value close to the 
conductance quantum and a minimum that becomes smaller for an increasing number 
of divacancies. In fact, conductance never crosses the critical value of one conductance 
quantum once a certain number of defects have been included. This result can be very 
clearly seen in Figure [IHl where the conductance never reaches values higher than 1 Go 
after N = 5 [7]. To clarify this issue, we have diagonalized the transmission matrix [30] 
to calculate the conductance for each one of the two conduction channels existing for 
energies around the Fermi level. According to this, only one of the two conduction 
channels is contributing to the total conductance once a number of 5-9 divacancies are 
presented in the nanotube. The origin of this channel suppression remains, however, 
unclear. 
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Figure 9. Conductance as a function of energy for the (5,5) nanotube and a mean 
distance between defects of d = 10.1 nm for a different number, N, of defects: N= 3, 
7, 10, 15 and 20 (from top to bottom). 



2 




number of defects 

Figure 10. (color online) Differential conductance at the Fermi level as a function of 
the number of defects, N, for d = 75.5 nm for the (10, 10) nanotube and for a selected 
number of the calculated random distributions of defects. 



The appropriately averaged resistance -according to the theory of disordered one- 
dimensional systems [29]- for a statistically significant number of random configurations 
of defects (up to 100 in the case of zero temperature, and up to 15 for the room 
temperature case, where the dispersion with respect to the mean value is much smaller) 
has been calculated for mean distances d between 10.1 nm and 75.5 nm. 

Figure [TT] shows the resistance, for the (10, 10) CNT, of a selected number of the 
calculated random distributions of defects at both zero and room temperature. Curves 
with the same color in (left) and (right) plots correspond to the same configuration 
of defects. At zero temperature, the resistance for each defect realization strongly 
fluctuates along the nanotube length, showing the characteristic behavior of disordered 
ID system. The average resistance, however, presents a clear exponential behavior as 
a function of the length, which is a typical feature of strongly localized systems. At 
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Figure 11. (color online) (10,10) CNT: Colored lines: ln{2R/Ro) as a function of 
the number of defects, N, for d = 75.5 nm and for a selected number of the calculated 
random distributions of defects at T = (left) and T = 300 K (right). Thick black 
lines: averaged resistance for the total number of defect realizations considered. Grey 
shadowed region: mean quadratic deviation with respect to the averaged resistance 
value. Ro is the inverse of the conductance quantum Gq. 



room temperature, resistance fluctuations disappear for increasing values of T, leading 
to very smooth curves when a critical temperature (scaling with 1/d) is reached, but 
preserving the exponential behavior. This fact conflrms the presence of the Anderson 
localization regime at room temperature for the defected nanotubes studied. This result 
has been conflrmed by experimental data [6,31]. 




5 10 15 20 
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Figure 12. (color onhne) (left) Averaged ln{2R/Ro) for different mean distances 
between defects d at room temperature for the (7, 7) nanotube. (right) Values for Lq, 
calculated by means of (p!Q|) . as a function of d for the same nanotube. 



Figure [I2] shows the averaged resistance for different mean distances at room 
temperature for the (7, 7) nanotube. For all defect densities studied we have found the 
strongly localized regime and the exponential increase of the resistance as a function of 
the tube length. We can then flt the average resistance for every d to the expression 



i?(L) = ^Ro e^/^° 



(10) 



where L is the length of the nanotube calculated as L = N- d and Lq is the localization 
length. In this way we can calculate the localization length Lq for every defect density. 
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These results for the (7,7) nanotube are summarized in Figure [T2l (right). A detailed 
analysis of the methodology and theoretical results obtained at different temperatures, 
as well as the experimental data supporting these conclusions, can be found in references 
[6, 7] for the (10, 10) and in [31] for the (5, 5). 

5. Conclusions 

In summary, we have analyzed, by means of a combination of ab initio simulations and 
linear-scaling Green's functions techniques, transport properties of armchair (5, 5), (7, 7) 
and (10, 10) carbon nanotubes with realistic disorder, namely mono- and di-vacancies. 
In a first step, first-principles simulations for each defect have been performed in order 
to determine the most stable reconstructions. The conductance profile for every defect 
has been then analyzed. Our results show that reconstructed monovacancies present an 
almost negligible drop of the nanotube conductance at the Fermi level compared to that 
of divacancies. As the lateral orientation of the divacancy is found to be about 2 eV more 
stable than the vertical one for all the nanotubes considered, only lateral divacancies 
have been included in our analysis of long tubes with different defect densities. Using 
an 0(N) recursive procedure based on standard Green's functions techniques we have 
calculated the localization length for several defect densities, predicting the existence of 
the Anderson localization regime for defected carbon nanotubes at room temperature. 
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